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Abstract 

We consider multiscale stochastic dynamical systems. In this arti- 
cle an intermediate reduced model is obtained for a slow-fast system 
with fast mode driven by white noise. First, the reduced stochastic 
system on exponentially attracting slow manifold reduced system is 
derived to errors of 0(e). Second, averaging derives an autonomous 
deterministic system up to errors of 0{y/e). Then an intermediate re- 
duced model, which is an autonomous deterministic system driven by 
white noise up to errors of 0{e), is derived using a martingale approach 
to account for fluctuations about the averaged system. This interme- 
diate reduced model has a simpler form than the reduced model on 
the stochastic slow manifold. These results not only connect averag- 
ing with the slow manifold, they also provide a martingale method for 
improving averaged models of stochastic systems. 



MSG: 34C15; 37H10; 60H10 

1 Introduction 



Complex dynamical systems in science and engineering often have multiple, 
disparate, time scales. Here we consider the common case of two widely 
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separated time scales. Recall that finite dimensional, deterministic systems 
are often written as the following slow-fast system of differential equations 

= Ax^ + f{x%f), (1) 

= \Bf + \g{x\f), (2) 

where e > is a small parameter, x^{t) G M" and y'^{t) G M*" are the state 
variables, A : M" ^ R" and B : MJ^ — > are bounded linear operators, 
and / : R" X R*^ ^ R" and ^ : R'^ X M"* ^ R'" are sufficiently smooth. 
For small e, the averaging principle suggests that a good approximation of 
the slow motion of the slow variables on long time intervals is obtained 
by averaging ([1]) over the distribution of the fast variables [Sj HI [HI e.g.]. 
The theory of invariant manifolds also supports, after the decay of exponen- 
tial transients, the reduction of the dynamics of the system ([I])-® to the 
dynamics of just x'^{t) [El HEl e.g.] which being lower dimensional is simpler 
to represent and analyse. 

For such deterministic systems, averaging is close to the slow manifold dy- 
namics. To derive a simplified system, the invariant manifold method seeks 
to construct a lower dimensional, smooth, attracting invariant set which gives 
a lower dimension reduced system by restricting the system to this manifold. 
If the gap in the spectra of A and B/e is large enough, then an exponen- 
tially attractive slow manifold in the form y = h{x) + 0{e) may be proved 
provided / and g are Lipschitz. Then, to errors typically of C(e), the lower 
dimensional, reduced system is 

X = Ax + f{x, h{x)). (3) 

On the other hand, in situations when the above slow manifold exists, the 
graph y = h{x) is the unique fixed point, exponentially stable, of the fast 
dynamics with frozen x G R". Then, by the averaging principle, an aver- 
aged system exists and indeed is the reduced system [9\. That is, in this 
deterministic case, if the system is autonomous, the leading approximation 
to the slow manifold reduced system coincides with the averaged system. 

However, in the stochastic dynamics explored in this article, the above close 
coincidence no longer holds. We explore the case when the fast variables y*^ 
are perturbed by a stochastic force modelled by white noise W; that is, 
replace (EJ by 

y^ = }:By^ + }:g^a:%y^) + ^W. (4) 

In this case the previously used fixed point of ([2]) instead becomes h'^{x, t), an 
ergodic stationary process of the stochastic dynamics (jlj) for frozen x G R". 
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Section [3] shows that averaging obtains, by replacing h{x) with h''{x, t) in ([3]), 
a random slow manifold reduced model. In this case the reduced system 
is a nonautonomous system with random stationary coefficient. Thus, the 
averaged system is a deterministic autonomous system with /(x, h{x))) in ([3]) 
replaced by 

1 r 

/(a;)= hm - / f{x,h^{x,t))dt. 

The ergodic property of h'^{x,t), Section [5l ensures that the effective non- 
linearity / is well defined. Since a random slow manifold, reduced system 
generically contains stochastic effects, necessarily the averaged system is dif- 
ferent. The two approximations contain different information of the stochas- 
tic dynamics of the original system. 

Analysing a stochastic system is considerably more difficult than a determin- 
istic autonomous system. Thus the main goal of this article is to derive a 
reduced model intermediate between the simple but deficient averaged model 
and the complexity of the random slow manifold model. The intermediate 
reduced model has errors 0(e) when compared to original system, whereas 
the averaged model has errors (9(e^/^), see section |5l A key step here is 
to estimate the fluctuation of the slow manifold reduced system from the 
averaged system, which is proved to be a Gaussian process by a martingale 
approach. 

Section [3] gives only the ffist order asymptotic approximation in e to the ran- 
dom invariant manifold. For an example. Section H] illustrates the results via 
a simple example model. Higher order terms in such an asymptotic approx- 
imation can be obtained by a normal form coordinate transformation |14j . 

SchmalfuB and Schneider [15] recently studied the slow manifold of a slow- 
fast random dynamical systems by a random graph transformation on a slow 
time scale 0{e). The random slow manifold is shown to be a family of random 
fixed points for the fast system (jl]) parameterized by x with e = 1 . But this 
approach does not give much information of the behavior of the slow system 
on the long time scales we aim for here. 
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2 Preliminaries 



Consider the following stochastic slow-fast system: 

= Ax' + f{x%y'), x'(0) = xoGM", (5) 
= -By^ + -g{x%y^) + ^W, /(O) = yo e M"^, (6) 

with matrix A : — > , matrix i? : — > M*", and nonlinear interaction 
functions / : x ^ M", ^ : x are C'^-smooth, k>2, with 

/(O, 0) = /x(0, y) = , (7(0, 0) = gy{x, 0) = . Here a is a nonzero constant 
and 1 > e > is a parameter. {W{t), t G M} is a two-sided M™- valued 
Wiener process. Moreover, / and g are Lipschitz with Lipschitz constant Lf 
and Lg respectively. Denote by | ■ |]Rn , | • and | ■ |m" + | ■ Ir™ the distance 
defined on M", M™' and M"+"^ respectively. For matrices A and B we assume 

Hi there are constants a > > (3 such that for any x G and y G M™" 

Ic^'^xIk" < e"*|a;|Kn , t<0, and |e*^?/|i8m < e'^*|2/|Km , t>0; 

(7) 

H2 the following spectral gap condition holds for any 1 > e > , 

for some A with P < eX < a . 

Under the above assumptions we give some basic results on {x''(t),y''{t)), 
the solutions of the slow-fast system (l5])-(l6l). First we have the following 
bounded estimation. 

Lemma 2.1. For any (xq, yo) ^ M""*"™, there is a unique solution {x'^{t), y'^{t)) 
to the slow-fast system such that for any T > 

(x^ y') G L\Q,C{0,T;W+"')) . 

Moreover there is a positive constant Ct such that 

E sup ?/'()!:)) |^„+™ < I (xo,?/o)|Rn+,n . 

0<t<T 

Proof. The existence and uniqueness result follows from a classical approach [1] . 
By the assumption on / we have 

^^k'(t)|Rn < \\A\\^,^\x\t)\ln + Lf\x'{t)\l^ , 



where ||v4||iru is the norm of A as a hnear operator on 
assumption on g and H2 we have by Ito's formula 



Also by the 



1 d 
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1 



a 



2dt 



Then by the Burkholder-Davis-Gundy inequality p!0] we have the bounded 
estimation result. This completes the proof. □ 

Now denote by yu'^ the probability measure defined by the distribution of the 
slow modes x*^ on space C(0,T; M") . Then we have the following tightness 
result. 

Theorem 2.2. The family of probability measure is tight in space 

C{0, T; M") for any T > . 

Proof By §^ 

x'[t)=xo+ [ Ax'{s)ds+ [ f{x'{s),y'{s))ds . 
Jo Jo 

For any < s < t we have 



E 



Ax^is) ds 



ElAx'is) 



ds 



1/2 



and 



E 



fix'{s),y'{s))ds 



< Vt 



E\fix'{s),f{s))\l„ds 



Then by the estimates of Lemma [2TT1 the assumptions on A and /, the family 
{x'^}e C L^(f2, C^/^(0, T; M")) is uniformly bounded with respect to e, which 
yields the tightness of {^^}e in C(0,T; M"). The proof is complete. □ 



3 Random slow manifold reduction 

In this section a random slow manifold reduced system is derived for the 
slow-fast system (E])-® for any small e > . We apply the theory of random 
invariant manifolds to study the asymptotic behavior of the system (E])-®- 
Theorem 13.91 proves that the long time behavior of (I5l)-(l6]) is described by 
the flow on a random invariant manifold which is exponentially stable. Then 
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by an asymptotic approximation, a random slow manifold reduced system is 
constructed. 

Recall some basic concepts in random dynamical systems (rdss) and random 
invariant manifolds. We work on the canonical probability space {Q, JF, P) 
with Q consisting of the sample paths of W{t). To be more precise, W is the 
identity on Q, with 



and P the Wiener measure [2]. 

Let Of : {Q, JF, P) —>■ (fi, JF, P) be a metric dynamical system (driven system), 
that is, 

• 9o= id, 

• OtOs = Ot+s for all s, t G M, 

• the map (t, uj) ^ OtOJ is measurable and ^^P = P for all t G M . 
On Vt the map 6t is the shift 



Definition 3.1. Let {X.,d) he a metric space with Borel a-algebra B, then a 
random dynamical system (rds) on {X,d) over 6t on (i7,jF, P) is a measur- 
able map 

If : R'^ X X X ^ X , (if:, cu, x) ^— s> ip{t, uj)x 
having the following cocycle property 



for t,T G R+, X E X and G i7 . 

A RDS (f is continuous or differentiable if (p{t,uj) : X ^ X is continuous or 
differentiable, respectively [2]. 

For a continuous random dynamical system ip on X given by Definition 13.11 
we need the following notions to describe its dynamical behavior [6j. 
Definition 3.2. A collection M = M{uj)i^^q of non-empty closed sets M{uj) C 
X , u E Q , is a random set if 



n = {we C([0, oo), M™) : w{0) = 0}, 



9tuj{-) = uj{- + t) -uj{t), teR 



G n . 



(9) 



(p{0,uj)x = x, (p(t,9rUj) o {p(^r,uj)x 



ip{t + r, uj)x 




is a real valued random variable for any x G X . 
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Definition 3.3. A random set B{u) is called a tempered absorbing set for a 
random dynamical system (p if for any bounded set K <Z X there exists txii-o) 
such that for all t > txi^j) 

ip{t,e.too,K) C B{uj), 

and for all e > 

lim e"^*(i(i3(6'_tcj)) = , for a.e. u e Q , 

t^oo 

where d{B) = sup^.gg (i(a;, 0), with E X , is the diameter ofB, ifOEB. 

Then we introduce the random invariant manifold [2, 7J. 

Definition 3.4. A random set M{uj) is called a positive invariant set for a 

random dynamical system (p(t,Lj,x) if 

!f{t,uj,M{uj)) C M{etuj), forty 0. 

If M{uj) = {xi + ip{xi,uj) : Xi G Xi} is the graph of a random Lipschitz 
mapping 

with X = Xi(B X2, then M{uj) is called a Lipschitz invariant manifold of ip. 

Duan et al. [7] give further details about random invariant manifold theory. 

For our purpose we introduce a stationary process r]'^{t) which solves the 
linear stochastic differential equation 

dr]'{t) = -Br]'{t)dt+ ^dW{t) . (10) 

Then by the assumption (Hi), ri^{t) is exponential mixing. Moreover, under 
the driven system we write the stationary process as ri^{9tuj). Denote by 
91 = 6t/^ the scaled version metric dynamical system 6t . Then, by Lemma 3.2 
proved by SchmalfuB et al. [I5], rj''{0tuj) has the same distribution of rj{6lu) 
with ri{t) = ri{9tuj) solving the following linear stochastic differential equation 

d'n{t) = Br]{t)dt + adW{t). (11) 

Now we consider the following slow-fast random differential system 

X' = AX' + F{X',Y% eiuj), X'{0)=Xo, (12) 

= -BY' + -G(X^ y^ etiu) , Y'{0) = Yo , (13) 
e e 
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with 



F{X^{t),Y^it),9tco) = /(X^(t),y^(t)+7^^(M), 
GiX%t),Y^{t),ep) = g{X%t),Y%t)+v%etUj)). 



(14) 
(15) 



By the assumption on / and g, solutions (X''(t), y (t)) to (fT2l)-(fT3l) define a 
continuous random dynamical system $'^(t,ct;) on M""*"™ defined by 

$^(t,cu)(Xo,Fo) = {X'{t,u),Y%t,u)). 

Then (p'^(t,uj) = ^^(t,Lj) + {0,ri^{9tu)) is a continuous random dynamical 
system defined by the slow-fast system 

Next we show that there is a random invariant set Ai{uj) for $'^(t,u;). We 
follow the method of Lyapunov-Perron for random dynamical systems [7]. 
For some P < X < a define the Banach space 



= {u:{-oo,0] 
with the norm 



pn+m . 



: u is continuous and supe \u{t)\]^n+m < oo} 



\u 



sup e ^^\u{t)\]^n 



te(-oo,o] 

For any given Xq G M" , define the map T{Xq, uj, ■) on by 
r(Xo,a;,(X^n)W = (e'^X, + j\'^'-'^^F{X\s),Y\s),eiuj)ds, 



,{t~s)B/e 



G{X'{s),Y'{s),eiuj)ds) . (16) 



Then, for any^Xo G M" , T{Xo,uj, 
(X^ Y'), (XS Y') G C;^ we have 



. Moreover, for each 



|T(Xo, iu, (XS Y^)) - T(Xo, cu, (X^ F^)) 



< sup 

t<0 



< 



< 



+ sup 

t<0 



[ e-^*e(*-^)^ [i^(^'(s), Y'{s), eiuj) - F{X'{s), Y'{s),ei 
Jo 

- / e~^*e(*-^)^/^ [G{X'{s), Y'{s), Olu) - G{X'{s), Y'{s), 9 

^ J —oo 



e 



,{f^-e\){t-s)/e 



UX',Y') - {X%Y' 



a — A 



|(x^y^) - {x%Y' 
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Then by Assumption H2, and the Banach fixed point theorem, for any Xq G 
R", has a unique fixed point (X", F") e with 

X'*{t) = e*^Xo+ / e^'-'^^F{X'*{s),Y'*{s),9luj)ds, 
= - / e^'-'^^^'G{X'*{s),Y'*{s),eiuj)ds. 



Now for each G , define a map cu) : M" M."^ by 

1 /"^ 

Xo ^ h'{Xo, uj) = - e-^'^'G{X'*{s), Y'*{s), 9liu) ds 
^ J —00 

which is Lipschitz continuous with Lipschitz constant 

J^h — T 



(17) 



{a-X){l-Lg [I /{a - A) + l/(eA -/?)]}' 

Then by a similar discussion to that of Duan et al. |7j we draw the following 
result. 

Theorem 3.5. Suppose assumptions Hi and H2. There exists a Lipschitz 
continuous random invariant manifold Ai'^{uj) for uj) which is Ai'^{uj) = 
{{X, h'{X,uj)) : X e M"}. 

Remark The random manifold A4^{uj) is independent of the choice of A 
which satisfies assumption H2. 

By the definition of (f'^{t,u!) we have the following corollary. 
Corollary 3.6. Suppose assumptions Hi and H2. The fast- slow system 
^ has a Lipschitz continuous random invariant manifold which is written 
as M'{uj) + {0,r]'{uj)) . 

Remark From ( JT7j) . by the transformation s/e s 

h'{X,uj) = / e-'^GiX'*{es),Y'*{es),eiiu)ds 
J —00 



in distribution) = f e-'^G{X'*{s),Y'*{s),esUj) ds {It 
J —00 



where 



X" (t) = e^'^X + e e^(*-")^F(X^* (s) , (s) , ^,cu) , t < 



Y'*{t) = I e^'-'^^G{X'*{s),Y'*{s),9sU;)ds, t<0 

00 



is the unique solution in C^^^ of the following system with a slow time scale 

X'it) = eAX'{t) + eF{X'{t),Y'{t), Otuj), (19) 
Y'{t) = BY'{t) + G{X'{t),Y'{t),etUj). (20) 

This can be deduced by the same approach to construct /i^(X, u). Then M.''{uj) 
is also an invariant manifold for system fll9l) - fl20l) driven by 6t. 

Now for small e > we give an asymptotic approximation for h'^{X,uj) for 
fixed uj & Q . For each X G M", consider 

Y{t) = BY{t) + G{X, Y{t), Otio) . (21) 

By the assumption on g and Hi, fl?T]) has a unique solution in Cj^ for any 
e > 0, denoted by Y{t) = h^{X, Otuj) , t < , with 



1-0 

'(X, uj)= e-'^G (X, /i°(X, Osuj), 9suj) ds . 

J — oo 



The existence of h^{X,uj) can be deduced by exactly the same approach to 
construct /i^(X, cj). Then we have the following approximation result. 
Lemma 3.7. Assume Hi and H2 hold. Then for almost all u & Q , 

\h'{X,uj)-h\X,uj)\R^ = 0{e) as e , 

which is uniformly for X in any bounded subset o/M". 

Proof By Remark El for any X G R'' , t < , 

|x^*(t)-xL _ 

rO 

< \e^'^-'^^'X - e-'^'Xl^n + eLf\{X'*,Y'*)\^,^ / e^-'^+'^^^^'-'Us . 

Then by the assumption Hi and H2, 

|X"(t) -X|^- = C(e) ase^O. 



Further 

|F-(t)-F(t)|^ 



Then also by assumption Hi and H2, for any X G MJ^ and G , 

\h'{X,u)-h%X,uj)W^ < \Y'*{t)~Y{t)\c-^ = 0{e), e^O. 
This completes the proof. □ 
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Remark Similarly, by Remark El we have 

1 

h%X,uj) = - e~^'^'G{X,Y'{s),9,/,uj)ds 

^ J — oo 

1 r° 

(in distribution) = - / e-^'^'G{X,Y'{s),eiuj) ds (22) 

^ J — oo 

where (X, F'(t)) G C^^ solves 

Y'{t) = ^BY'{t) + (X, Y'{t), eiuo) . (23) 

Then M.^{uj) = {{X , {X , u)) : X G M"} is an invariant manifold for sys- 
tems (EH) and ( 1231) driven by 6t and respectively. 

Given a random invariant manifold, in the following, if the spectral gap is 
large enough, any solution of (|T^ - (fT5|) is proved to converge exponentially 
quickly to a flow on the manifold Ai'^{uj) as t ^ oo . Then the {n + m)- 
dimensional slow-fast system (l5])-(l6l) is reduced to an n-dimensional system. 
First we give the following definition for a random dynamical system $. 
Definition 3.8 (Almost sure asymptotic completeness). Let Ai{uj) be an in- 
variant manifold for random dynamical system $(t, uj). The invariant mani- 
fold M. is called almost surely asymptotically complete if for every z G M^^™, 
there exists z' G M.{uj) such that 

\^{t,uj)z -^{t,uj)z'\ < D{uj)\z - z'\e-^\ t > 0, 

for almost all uj E Vt, where k is some positive constant and D is a positive 
random variable. 

Following the approach of Wang and Duan [19] we have the following result. 
Theorem 3.9. Assume Hi and H2. If there is 5 > Q such that 

ea + eLf + e5Lf + 13 + Lg + <Q , (24) 

then the random manifold M.{uj) is almost surely asymptotic complete. That 
is, for any solution (X^(t, cj), y^(t, cu) of ( ITgl) -( IT^) . there exists an orbit 
(X'^(t, cu), y^(t, u;)) on the manifold M.{9tuj) which is governed by the fol- 
lowing n- dimensional different equation 

X' = AX' + F{X%Y',eiu:), Y' = h'iX^OtUj), (25) 

such that for any t > and almost sure uj E Q 

I (X^(t, 00), Y'{t, 00)) - (X^(t, 00), Y'{t, a;)) 
< D\{Xo,Yo) - (X^(0),F^(0))|M.+™e-^*/^ 

with 7 = — /5 — Lg — S^^Lg > and some deterministic constant D > . 
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Remark An important problem is to predict the initial value (X'^(O), i^'^(O)) , 
which one can deduce from a stochastic normal form [T3] or related dynamical 
considerations [13]. 

By the assumption on /, from the reduced system (l25l) and Lemma [3?7l we 
have 

X' = AX' + F{x\h'{x\eiu),eiuj) 

= AX' + F{X\ h\X% eiu) + C(e), 61^) 
= AX' + F{X% eiuj), e'^uj) + 0{e). 

Then we draw the following slow manifold reduced result. 
Corollary 3.10. The slow manifold reduced system for the slow-fast sys- 
tem up to errors ofO{e), is 

x'it) = Ax'it) + f{x'it), h\x\t), etuj) + viOtuj)) . (26) 



Remark On the other hand, by assumption H2, Lg (3 < which yields 
that for any fixed X G M", system ( |T3|) has a unique stationary solution 
which is, by the definition, co"). Moreover, the stationary solution is 

exponentially stable for almost all u; G . So the slow manifold reduced 
system fpUj) is obtained by replacing y' by y'^^^u) = h^{x,u!) + rj'{uj) the 
stationary solution of for fixed x. This is the same as the deterministic 
autonomous case. 



4 Application to a simple example model 

Consider the following nonlinear slow-fast stochastic system with {x, y) G 



dx = f{x,y)dt, (27) 

dy = -[-y + g{x,y)]dt+^dW{t), (28) 

where / and g are smooth with the following specific values close to the 
origin, and far from the origin 



f{x,y) 



-xy , + < 1/16 
, + > 1/8 , 
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and 

jx^-2y\ x2 + |/2<i/i6, 
10, + > 1/8 . 

In a slow time scale we rewrite the above system as 

dx = ef{x,y)dt (29) 
dy = [-y + g{x,y)]dt + adW{t). (30) 

The nonlinearities / and g satisfy assumptions Hi and H2. Then by The- 
orem [33], the example system (l27j) -( |28ll has an invariant manifold Ai{u!) = 
{{x, h'^{x,u) + ri'^{uj)) , X G M} and the reduced system on this manifold is 

dx = f [x, h'{x, Otuj) + r]'{etuj)\ dt . (31) 

Further by the definition of / and g and ( l22l) , h'^ has the following asymptotic 
expression near x = 



r"(0) = hQ{x,uj) + 0{e) 

e'/'[x^ - 2{ho{x, eioj) + r]{eiuj)f]ds + 0{e) 



00 



,x^] 



- f e'/'[x^ -2{x^ + r]{eiu)f]ds + 0{e, 
^ J —00 


e'''[x^ - 2{x^ + r]{9lu)Y]ds + C(e, x^) 

/ e'/'[r]{9luj)]'^ds / e'''7]{eiuj)ds - 2x^ 

J —00 ^ J —00 



2 2 

X 

e 



—00 
3^ 



+ 0(e,x 

Then an asymptotic expression of the reduced dynamics of the example sys- 
tem (EH) is 



x = -x^- X7]\etuj) + — / e'/'[r]\e,+tuj)fds + / e'''7]\et+sUj)ds 

+ 0(e,x^). (32) 

For small e the random integrals in ( l32l) can be replaced by some stochastic 
term on a long time scale as follows [5, e.g.]. 
Lemma 4.1. Fore small, in a mean square sense 

-^r]^{t) -^adWiit), f e'/'r]%e,^too)ds ^adW^it), 
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and 



with some mutually independent standard scalar Brownian motions Wi{t), 
W2it) and Wait). 

Proof. This follows by a martingale approach. Let 

1 



Ve Jo 

then z"^ solves 

z' = ^r^\t), z^(0) = 0. 

Then by martingale results [211 EDI e.g.], we have the limit of as e — > 0, 
solves 

dz = a dWi , ^(0) = , 

for some scalar Brownian motion Wi (t) . The others can be obtained similarly. 

□ 

Thus for small e, the reduced system f lHTl) is 

dx ={-x^ + a'^x)dt - y/eaxdWi{t) + V2ea^xdW3{t) + A^eax^dW2{t) 
+ 0{e,x^). (33) 

Observe that the three noise terms are all 1/2-order in e; further, the three 
noise terms can be replaced by one noise term. This will be treated by 
averaging in next section. 

Note that the reduced model ( l33l) is same as that deduced by the stochastic 
normal form [14j of the slow-fast system ( l27|) -(l28l). 



5 Averaging approximation for small e 

The random invariant manifold reduction gives a lower dimensional system 
which is a nonautonomous random system. This section, for small e, derives a 
simpler reduced system which is an autonomous system perturbed by a small 
stochastic term. Moreover, the approximating errors prove to be of 0{t). 
First, as e is small in the slow-fast system (I5])-(I6]), y*^ is highly fluctuating in 
time, and so an averaging method derives the averaged equation ( l34l) which 



14 



approximates x*^ up to errors of 0{y/e). Second, by a martingale approach, 
we derive an intermediate reduced equation (jH]) which approximates x"^ to 
errors of 0{e). 

We just need to consider the reduced system (l26l) . First define the following 
averaged equation 

x{t)=Ax{t) + f{x{t)), x{0) = Xo, (34) 

with the averaged 

/»=E/(x,r'"(^)) 

which is independent of e. By the Lipschitz property of h'^{x,uj) and the 
assumption on /, the average / is Lipschitz; denote the Lipschitz constant 
by Lj. Moreover, by the exponential mixing of y'^'^ we have for any solu- 
tion y^'^{t) of (El) with initial value yo G and any fixed a; G M" 

nf^^it) - r'^'(t)lR"^ < Ebo - r'^'(0)re^*/^ . (35) 

Here y'^'^{t) = y'^'^{6tuj) . Then we prove the following result 

Theorem 5.1. Assume Hi and H2. Given any T > , for any Xq G M", as 

e — s> the solution of ^ converges in probability in C(0,T;R") to x which 
solves [34\ )- Moreover, the rate of convergence is 1/2, that is, 

sup E\x'{t) - x(t)|Mn < CtV^ 

0<t<T 

for some positive constant Ct- 

Proof. For any T > , we partition [0, T] into subintervals of length 5 = 
y/e. Introduce a new process x''(t) satisfying for t G [kS, {k + 1)5) , k = 

0,1,... ,[r/5], 

F(t) = e(*-'=^)V(M) + [ e^'-'^^f [x'{k5),y'^''^^^^\s)]ds. 

Jkd 

Then for any t G [kS, {k + 1)6) 

E\x%t) - X%t)\^r. 

< E [ e(*-^)^ 1/ [x^ik6),r'^'^''\s)] - f [x^is),r'''^^\s)] ds 

JkS 

< E [ e^^~'^^Lf[\x'{k6) - x'{s)kn + \ f'-'('^'\s) - r'^'^^'\s)W^] ds 

JkS 

= E [ e^'-'^^Lf[\x'{k5) - x'is)]^^ + \h'{x'ik5),9sio) - h'ix'{s),9suj)\^^] ds 

JkS 

15 



JkS 

< CtS 

for some positive constant Ct- 

On the other hand by and noticing the Lipschitz property of / we have 

E|x^(t) -x(t)|Mn 

Jo 



f[x'{[s/6\)]-f[x^{s)]L ds 



3(*-^)^|/[x^(.)]-/[x(.)]L„ ds 



e+ / E|x''(s) — x(s)|iRn 



+ E 
+ E 
< Ct 

for some positive constant Ct- Then we have 

E\x'{t) - x{t)\^n < E\x'{t) - F(t)|Mn + E\x'{t) - x(t)|i 
And by the Gronwall lemma we have for any t G [0, T], 

E\x'{t)-x{t)\Kn < CtV^. 
This completes the proof. 



□ 



Although the averaged equation is simple — being a deterministic autonomous 
system — the above result shows the approximation error is of 0{^/e). Now 
in the following we give a refined approximation of x*^ to errors of 0{e) for 
small e > . For this define 

i/(a:,t/e) = /(x,r'^(t))-/(x). 

Then rewrite (l26l) as 

x'{t) = Ax'{t) + H{x'{t),t/e) + f{x'{t)). (36) 

We follow a martingale approach to prove that H{x''{s), s/ e) ds can be 
approximated by a Gaussian process. Define 



H{x,t) 



E[H{x,t + s)\Tt] ds 
16 



E [H{x, s) I Tt] ds (37) 



and 



1 /■* 

= ^/^[H{x'{t),t/e) - H{xo,0)] + ^ H{x'{s),s/€) ds (38) 

Jo 

-yfe I H^{x'{s),s/e) [Ax'{s) + H{x'{s), s/e) + f{x'{s))]* ds 
Jo 

where * denotes the transpose of a vector. Then we have the following result. 
Lemma 5.2. For each e > , {M^ : <t <T} is a martingale with respect 
to {J-'t/e : ^ > 0}, and the quadratic covariance is 

{M')t = e[H{x'{t),t/e)Y -e[H{xo,0)Y -2^ f H{x'{s), s/e) dM'^ 

Jo 

-2e [ H{x'{s), s/e)H^{x'{s), s/e) [Ax'{s) + H{x'{s), s/e) + f{x'{s))] * ds 
Jo 



+ 2 / H{x'{s),s/e)H*{x'{s),s/e)ds. 



(39) 



Proof. Let [^1,^2] be an arbitrary subinterval of [0,T], A = {soj Si, . . . , s^v} 
a partition of interval [^1,^2] and h = maxi<j<jv{si — Si-i} • Then 

E lH{x'{t2),t2/e) - H{x'{ti),ti/e) + ^ j[ ' H{x'{s), s/e) ds 



t2 



H^{x'{s), s/e) [Ax'{s) + H{x'{s), s/e) + f{x'{s))] ds 
= ^{Y.r [Hx{x\s), s^/e) - H,{x\s), s/e)] 

y i=i Jsi-i 

X [Ax'{s) + Hix%s),s/e) + f{x'{s))] 

/ ^ [H{x'{si_r),s/e)-H{x'{s),s/e)] ds 

By the assumption on / and the definition of H we have 

^[\Hx{x''{s),s/e) - H^{x^{s),Si/e)\^„^Tti} <Lfh, se [si-i,Si\, 



ti/t 
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and 

\Hix'{s,.,),s/e)-H{x'{s),s/e)\^„ 

+ |/(x^(3,_i),r'^^W(s))-/(x^(.),r'^^W(.))|^„ 

< 2Lf{l + Lh) \x\si.i) - x\s)\^^ . 

Then by the estimate in Lemma ET] and the Jensen inequahty for conditional 
expectation, passing to the hmit /i ^ yields the first statement of the 
lemma. 



By the Ito formula 



Then 



d[Mlf = 2MldMl + d{M') 



t ■ 



{M')t = [M,^]^ -2 [ Ml dMl . 
Jo 



Substituting the expression for M* into the above equation and integrating 
by parts yields flH^ . The proof is complete. □ 



Now by the definition of , we further rewrite fl26p as 

x'{t)=xo+ [ [Ax'{s) + f{x%s))]ds + y/^M^ + eR'{t), (40) 
Jo 



where 

R'{t) = -H{x%t),t/e) + H{xo,0) 



+ [ H^{x'{s),s/e) [Ax'is) + H{x'{s),s/e) + f{x'{s))Y ds . 
Jo 

By the estimate of Lemma I^TTl the assumption on /, and the definition of H, 

0. (41) 



limE 



sup y/eR'{t) 

0<t<T 



Now define the process 

1 
7^ 



Ml = ^{ x\t) -xo- I [Ax'is) + f{x'{s))] ds \ . (42) 
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By the exponential mixing property of y'^'^{t), the assumption on /, LemmaEU 
and (HOD, {Ml}o<e<i is tight in space C(0,T;M"). Let P be a hmit point of 
the family of probability measures of {£(A1^)}o<e<i . And denote by a 
C(0, T; R")-valued random variable with distribution P. Let be a contin- 
uous bounded function on C(0,T;M"). Set *'(s) = ^(x'(s)), then by 
and fHTj) and Lemma [5.21 

E [(Ml) - M\)) ^\s)\ = E [{V~eR'{t) - V^R^s)) ^>%s)\ ^ , e ^ , 

which yields that the process is a P- martingale with respect to the Borel 
cr- filter of C(0,T; M"). Further by Lemma [5.21 we have 



E 


sup 


(MOt - 2 / 




.0<t<T 


JO 



0, 



0. 



Denote by the matrix 2E [//"(x, t)i7*(x, t)] , x G M", which is indepen- 
dent of t. By the definition of H, H and the exponential mixing property 
of 1/*^'^, we derive the following representation for S(a;): 



2E 
2 



E[H{x,s)\rt] ds H*{x,t) 
E[H{x,s + t)H*{x,t)]ds 
E[H{x,s)H*{x,0)]ds. 



(43) 







Notice that H{-,t) and H{-,t) are Lipschitz uniformly with respect to t. So 
we follow the approach used by Pardoux et al. \TT, Proposition 8]: 



e/sup [ [H{x'{s),s/e)H*{x'{s),s/e) -^{x'{s))]ds 

lO<t<T Jo 



0, e^O. 



Then we pass limit e ^ in ( H2l) and arrive at the following result. 
Lemma 5.3. The process 



M = Mm Ml = lim ^ <^ x'(t) - xq 



[Ax'{s) + f{x'{s))] ds 



defined on the probability space (C(0, T; M"), ;B(C(0, T; R")), P) is a square 
integrable martingale with the associated quadratic variation process given 
by S(x(t)) which is defined by (4^- 
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Proof. It is just necessary to prove the uniqueness of the martingale prob- 
lem. By the martingale representation theorem [TT], without changing the 
distribution of J^l and A^°, one can enlarge the original probability space 
to {fl, JF, jFj, P) and there is a ]R"-valued Wiener process W(t) such that 
satisfies 

dz{t) = a{x{t))dW{t), z{0) = , 

where a{x) G M"^" such that a{x)d-*{x) = x{t) solves ([34]). By the def- 

inition of E and H, a is Lipschitz which yields the existence and uniqueness 
of P . This completes the proof. □ 

Now in order to derive the intermediate reduced system we introduce z'^ 
which solves 

dz'{t) = a{x'{t))dW{t), z'{0) = . 

Then by the Burkholder-Davis-Gundy inequality [10] and Theorem 15. we 
have for any T > that 

E sup \z'{t) -z{t)\l„ < L^Cxe 

0<t<T 

for some positive constant Ct- Here E is the expectation operator defined on 
the enlarged probability space with respect to probability measure P. Then 
by the definition of Ai^, without changing the distribution of a;*^, over the 
enlarged probability space {Q, JF, P) we have 

x'{t)=XQ+ f [Ax'{s) + f{x'{s))]ds + ^ [ a{x'{s))dW{s) + 0{e). 
Jo Jo 

By discarding terms of 0{e) we derive the following reduced approximation 
system 

F(t)=xo+/ [Ax'{s) + f{x'{s))]ds + Ve [ a{x'{s))dW{s), 
Jo Jo 

or, equivalently, the following differential form 

dx'{t) = [Ax'{t) + f{x'{t))] dt + y^a{x'{t)) dW{t), x'{0) = xq . (44) 
Then we draw the following result 

Theorem 5.4. Assume Hi and H2. Over the enlarged probability space 
(f2,jF, P), without changing the distribution, for any T > , a;o G M", 
the solution of x", is approximated by x^ , the solution of (^^P, in the 
space C(0,T;M") up to errors ofO{e), that is almost surely 

- x^\c{o,T;R") = 0{e), e^O. 
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Remark By the above result, the system — which is a deterministic 
equation with small stochastic perturbation — is an intermediate reduced sys- 
tem for the original slow-fast system ©-(E]). Moreover, the system (jHj) 
approximates original system (I5])-(I6]) up to errors of C(e). 

For an example, we now return to the simple model of Section HI By the 
definition of / and g in Section HI for any fixed x G M , fl^Hj) has a unique 
stationary solution y'^'^ which is exponential mixing. Then near x = we 
have the asymptotic expansion 

And the averaged equation near a; = has the following asymptotic form 

x = f{x) = -x^ + a^x + 0{x\^/e). (45) 
To obtain an intermediate reduced system we calculate, near a; = , 

POO 

2E / [ - xy'^%t) + xEy'^%t)] [ - xy''%0) + xEj/^'^(0)] dt 
Jo 

/■oo 

= 2x' / [E(r'^(t)r'^(0))-E(r'^)']rft 
= xV2-12(tV + 20ctV + O(x3). 
Then we obtain the following reduced system 

dx = [-x^ + a^x] dt + (TXv^Vl - 12x2 + 20^2 dW{t). 

6 Conclusion 

Averaging of systems with noise in only the fast variables leads to a deter- 
ministic reduced model with errors 0{y/e), Section [5l The stochastic slow 
manifold. Section [3l shows that generally there should be noise in the reduced 
model, as seen in the simple example of Section HI Consequently, Section [5], 
shows that a martingale argument establishes the leading influence of the fast 
noise on the averaged system. When added to the averaged system the resul- 
tant reduced model typically has the smaller error 0{e). This work not only 
shows how to improve averaged models of stochastic systems, but strengthens 
the connections between the methodologies of averaging and slow manifold 
reduction. 
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